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ABSTRACT 

CN i We present a model for the spindown of young radio pulsars in which the neutron 



o 



star loses rotational energy not only by emitting magnetic dipole radiation but also 



^> ' by torquing a surrounding disk produced by supernova fallback. The braking index 

^ ■ predicted in our model is in general less than n = 3 (the value for pure dipole magnetic 

' radiation), in agreement with the reported values of n < 3 for five young radio pulsars. 

With an accuracy of 30% or better, our model reproduces the age, braking index and 
third frequency derivative of the Crab pulsar for a disk mass inflow rate in the range 
3 X 10^6 - lO^"^ g s-i. 

Q . subject headings: accretion, accretion disks - supernovae: general - pulsars: general - stars: 

|H I neutron 

' 1. Introduction 

5t , The magnetic dipole model for radio pulsars has proven extremely successful in explaining 

the observational properties of these objects (Pacini 1967; 1968; Gunn & Ostriker 1969; Gold 
1968; Goldreich &; Julian 1969; see, e.g., Shapiro & Teukolsky [1983] for an overview; Manchester 
& Taylor [1977] for a detailed review). A longstanding problem of the theory, however, has been 
that the predicted braking index (see §2.1 for definition) is n = 3, while measurements of braking 
indices in five young radio pulsars are all less than 3 (see §3 for the most recent measured values). 
Various effects have been proposed to explain this discrepancy, such as magnetic axis wandering 
(Macy 1974), non-standard vacuum dipole theory (Melatos 1997), and non-dipolar field structure 
or field strength secular evolution (Blandford, Applegate & Hernquist 1983; Blandford & Romani 
1988). 
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In this Letter, we explore the consequences of a recent proposal by Chatterjee, Hernquist & 
Narayan (2000) according to which young radio pulsars could be surrounded by remnant disks 
produced by supernova fallback (see also Michel & Dessler 1981, 1983; Michel 1988; Yusifov et al. 
1995; Alpar 1999, 2000). The disk may indeed be able to contribute to the neutron star spindown 
even in the radio pulsar phase, provided that matter does not fall through the magnetosphere 
inside the light cylinder radius. We note that Marsden, Lingenfelter & Rothschild (2001a, b) have 
recently considered some of the consequences of the presence of such a disk around a radio pulsar, 
with a special emphasis on explaining discrepancies between a pulsar's true age and its timing age 
(see also Gvaramadze [2001] for the case of an interaction with dense circumstellar material). The 
outline of this Letter is as follows. In §2, we derive a theory for pulsar spindown in the presence of 
a fallback disk and the associated propeller torques. We apply the theory to the Crab pulsar and 
other young radio pulsars in §3. 



2. Pulsar Spindown Model 

We certainly lack a proper understanding of the magnetospheric interaction between a radio 
pulsar and a surrounding disk. However, simple arguments can be used to derive a reasonable rate 
at which the pulsar is spun down by this interaction. Implicit to our model is the assumption that 
the disk can survive the high-energy environment created by the intense relativistic pulsar wind. 
For the pulsar mechanism to operate, we require the disk to be located beyond the light cylinder 
radius, at i?/c = c/Q (where is the neutron star angular speed and c is the speed of light). On 
the other hand, the magnetic field lines are open beyond Ric and the field structure is no longer 
dipolar. This is the main motivation for our assumption that the disk extends down to Ric- At 
Rlc, the inflowing material encounters a dipolar field structure rotating at nearly the speed of 
light: an efficient propeller effect (Illarionov & Sunyaev 1975) may be expected to operate there, 
spinning down the neutron star and preventing mass inside the pulsar magnetosphere at the same 
time.0 

Following Menou et al. (1999; see also Daumerie 1996), we adopt an efficient propeller torque 

j = m = -2MRiQ, (1) 

where / is the moment of inertia of the neutron star, a dot denotes a time derivative and M 
is the disk mass inflow rate at the magnetospheric radius, Rm- This expression for the torque 
assumes that the material flung away by the propeller effect has been accelerated to an angular 
speed corresponding to that of the star. In the pulsar case, we effectively take Rm = Ric and the 
material flung away only approaches the speed of the light. Under these simple assumptions, the 



■^We note that, if the disk were truncated beyond Ric (e.g. because of a large energy density right outside the 
Ught cyUnder), the torque on the neutron star would likely be largely reduced or entirely suppressed. 
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torque for this mildly relativistic propeller is given by 

j = -2Mc'^n-^, (2) 

which is essentially what one expects from dimensional arguments assuming the torque acts at 
A more efficient, relativistic propeller is energetically possible; it would correspond to a reduced 
value of the disk mass inflow rate for a given value of the torque. 

The pulsar loses rotational energy at a rate: 

E = mh = -pn^ - 7, (3) 

where the first term on the RHS corresponds to magnetic dipole radiation losses and the second 
term to propeller losses. We have 

7 = 2Mc' ^ 2 X 10^^ erg s-^ (yW?F^) , (4) 
^^sin^ aR%s 

P = ^^3 ' (5) 

where Rns is the neutron star radius and Bp is the magnetic field strength at the magnetic pole, 
inclined at an angle a from the rotation axis (see, e.g., Shapiro & Teukolsky 1983). 

In what follows, we assume that Eq. (^) is valid at all times. We note that this assumption 
is not entirely consistent with the evolutionary scenario of Chatterjee et al. (2000), in which a 
neutron star decouples completely from its surrounding fallback disk when it enters the radio 
pulsar phase. We expect this discrepancy to be minor given that magnetic dipole radiation still 
dominates the losses in the present model. More importantly, a young neutron star is not initially 
in the radio pulsar phase in the Chatterjee et al. model, but in a "dim propeller" phase (during 
which the magnetospheric radius has not yet reached the light cylinder radius). This phase can 
last up to several hundreds of years in the Chatterjee et al. model and it is not captured by our 
Eq. dll). This discrepancy suggests that the results derived below for the integrated neutron star 
spin history should be interpreted with caution (instantaneous quantities such as the braking 
index are not directly affected, however). 



and 



2.1. Steady mass inflow rate 

From Eq. (^), we find that the braking index is 

n=—. — = 3—- 3 , (6) 

if 7 (i.e. the disk mass infiow rate) is constant with time. The breaking index is therefore less than 
3, the value corresponding to pure magnetic dipole radiation losses (recovered in the limit 7 — > 0). 
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Following the usual procedure, Eq. @ can be separated and integrated to constrain the 
neutron star spin history. We rewrite Eq. (^) as 



^_27 

1 + ^04 / 



dt, 



(7) 



which can be integrated (assuming 7 = constant) by noting that the LHS is of the form 
dy/{\ + y^) = (i(arctan(?/)). This yields, at any time t after the system's birth: 



2(7/3)1/2 



arctan 



7 



arctan 



7 



(8) 



where rt{t) is the angular rotation speed at time t and is the angular rotation speed at birth. 
In the limit of very large the left arctan is simply equal to 7r/2. There is no simple relation 
between the pulsar timing age and true age in this solution, contrary to the case of pure magnetic 
dipole radiation losses (Shapiro & Teukolsky 1983; Gunn & Ostriker 1969). 



2.2. Time-dependent mass inflow rate 

We generalize here the above calculations to the case of a mass inflow rate decreasing with 
time, as expected in the fallback disk scenario (Chatterjee et al. 2000; Menou, Perna & Hernquist 
2001). Allowing the mass inflow rate to vary with time, the braking index becomes: 

1 _ lifl 

n = 2,-- (9) 

4 ^ 47'' 

where both 7 and (l are negative (assuming a decreasing mass inflow rate).^ The quantity —il/tl 
is twice the usual pulsar timing age (= P/2P), while the quantity —07/7 is actually the true age 
of the system if, as expected, the mass inflow rate obeys a power law evolution in time of the form 

7(t) oc 

While the braking index and the second deceleration parameter (defined in §3.1) are 
instantaneous quantities which can easily be computed for any power-law solution (using, if 
known, the values of the true age and the timing age as explained above), the spin history cannot 
be computed analytically for an arbitrary power law index a, but only in some special cases. 

For simplicity, here we restrict our discussion to solutions for the neutron star spin history 
with a time-dependent mass inflow rate satisfying ^{t) = 7o(i/^o)~^- Let w = l/fi^f, a = 2/3// 
and h = 2^01^/1, then the differential equation (^) becomes 

dw ,9 , , 

w + t — = a + bw^ (10) 



^We note that this expression for the braking index can potentially yield values of n > 3. 
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This can be separated and integrated: 

dw 



I 



Imp' — w + 



If Aah < 1 then the integration yields 

2bw - 1 - c 



c 



2bw - 1 + c 



Int = constant (11) 



— Int = constant (12) 



where c = \/l — 4a6 . The constant is fixed by the value of the variable w at some initial time ti. 
The solution for CI can then be easily derived: 

^2 ^ 26 {2lm',- l + c)-i2b,c,- l-c){l/l,:y' 

t {l + c){2bwi-l + c)-{l-c){2bwi-l-c){t/tiy ' ^ ' 

At later times, t ^ U, and the solution can be well approximated by 

n^^^. (14) 

2at ^ ^ 

Evaluating the above equation at the present time to yields 

to = -j^ ^. (15) 

Note that the factor 1/(2(50,^) is the pulsar timing age (in the case of pure magnetic dipole 
radiation losses). Although this solution corresponds to rather simple evolution laws for CI and 7, 
it is useful as an illustration of the effect of allowing for a time-varying mass inflow rate in the 
equations (as compared to the steady case considered in §2.1). Given that a fallback disk should 
have a mass inflow rate decreasing approximately as (Cannizzo et al. 1990; Chatterjee et 
al. 2000; Menou et al. 2001), the steady and simple time- varying cases considered here may be 
expected to bracket the real case (7 oc t~^'^), for which a numerical integration is required in full 
generality. 



3. Applications 

3.1. The Crab Pulsar 

For the Crab pulsar (PSR B053H-21), we adopt I « 10^^ g cm'^ and /3 = 3.85 x 10^9 erg s^. 
These values are derived from measurements of O, the estimated value of — £^ 5 x 10^^ erg s~^ 
and the assumption of pure magnetic dipole radiation losses (see Shapiro Sz Teukolsky 1983; 
Taylor, Manchester & Lyne 1993). The propeller losses are comparatively small in the applications 
considered below and for the Crab pulsar: 

^.2J J' X\ (16) 
7 Il0i6gs-M ^ ^ 
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We recall here that, in 1972, the true age of the Crab pulsar was 918 yrs, while its timing age was 
measured to be 1243 yrs (which is half the characteristic age as defined by Shapiro & Teukolsky 
1983). The braking index of the Crab pulsar is n = 2.509 (e.g. Lyne et al. 1988). 

In the case of steady mass inflow in the fallback disk outside the light cylinder, we find using 
Eq. (^) that the braking index is n = 2.85, 2.57 and 1.86 for the mass inflow rates M = 10^^, 
3 X 10^^ and 10^^ g s~^, respectively. For the same mass inflow rates, Eq. (^) predicts a system age 
of to = 1130, 1101 and 1018 yrs, respectively (in the limit of very large initial angular speed Oj). 

In the case of a mass inflow rate decreasing with time as we find using Eq. (|9|) that the 

braking index is n = 2.97, 2.93 and 2.78 for the (current) mass inflow rates M = 10^^, 3 x 10^^ and 
10^^ g s^^, respectively. The age cannot be simply computed with the power law index a = 1.2. 
However, it will be bracketed by the values obtained for steady mass inflow above and the a = 2 
case, for which Eq. ( |T5| ) holds. For the same (current) mass inflow rates as above, Eq. ( [Tsl) predicts 
a system age of to = 1099, 1020 and 817 yrs, respectively. 

In each of the cases above (steady mass inflow and 7 oc t~^), it is possible to find a mass 
inflow rate for which both the Crab pulsar true age and braking index are predicted to within 20% 
or better. Consequently, it appears that a fallback disk with a current mass inflow rate in the 
range 3 x 10^^ — 10^'^ g s~^ can account for the properties of the Crab pulsar with an accuracy of 
20% or so. 

The model can be further tested with the third time derivative of the Crab pulsar frequency 
(Lyne et al. 1988; Blandford &: Romani 1988). It is useful to define a second deceleration 
parameter: 

002 / 1 \^ i + iOi-iO 

^- 03 1 , + , , /304 > 

which has the correct limit p = 15 for 7^0 (pure magnetic dipole radiation). The expression 
for steady mass inflow is found by forcing 7 = 7 = in the above equation. Note again that 
a(a + 1)7/7 is the square of the pulsar true age for a power law evolution of the disk mass inflow 
rate (7(t) oc t^"), while —Vt/Vl is twice the pulsar timing age. The measured value of p for the 
Crab is « 10.2 (Lyne et al. 1988; Blandford & Romani 1988; Melatos 1997). We find, for the 
case of steady mass inflow, that the equation above predicts p = 12.6 for M = 3 x 10^^ g s~^ and 
p = 9.46 for M = 10^'' g s~^. For a time dependent evolution with a = 1.2, we find p = 13.6 
and p = 11.96 for the same mass inflow rates, respectively. The second deceleration parameter is 
therefore correctly predicted by the model (with an accuracy of 30% or so) for fallback disk mass 
inflow rates in the range 3 x 10^^ — 10^"^ g s~^. 
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3.2. Other Young Radio Pulsars 

Braking indices have been measured for four other young radio pulsars. The measured values 
are n = 2.2 ±0.1, 1.4 ±0.2, 2.91 ± 0.05 and 2.837 ± 0.001 for PSR B0540-69, PSR B0833-45 (Vela), 
PSR J1119-6127 and PSR B1509-58, respectively (Lyne, Pritchard & Graham-Smith 1993; Deeter, 
Nagase & Boynton 1999; Lyne et al. 1996; Kaspi et al. 1994; Camilo et al. 2001). 

It is possible to calculate, within the framework of our pulsar spindown model, the disk mass 
inflow rates required in each system to account for the measured breaking indices. In the absence 
of known true ages for these four pulsars, we calculate theoretical braking indices in two limits: (1) 
the steady mass inflow case (Eq. ^ and (2) the time-dependent case (Eq. |^ with the additional 
simplifying assumption that the pulsar timing age is a correct measure of the pulsar true age (the 
validity of this assumption, appropriate for the magnetic dipole model, was recently challenged by 
Gaensler &; Frail 2000). The pulsar properties are taken from Taylor et al. (1993) and Camilo et 
al. (2001). 

We find that the pulsar spindown model can account for the braking indices of PSR 
B0540-69, PSR J1119-6127 and PSR B1509-58 for disk mass inflow rates ~ 3 x 10^^, ^ lO^"^ and 
^ 10^^ g s^^, respectively. In each case, the propeller losses represent a reasonably small fraction 
of the total pulsar rotational losses. The very small value of n ~ 1.4 for the Vela pulsar cannot be 
accounted for by a model in which propeller losses are small. We note that the sub-Eddington 
mass inflow rates requhed for PSR B0540-69, PSR J1119-6127 and PSR B1509-58 (as well as 
the Crab) appear consistent with the values expected at late stages (before the fallback disk 
becomes neutral) in the evolutionary scenario described by Menou et al. (2001). One expects 
the magnetospheric interaction to be largely reduced or even suppressed when the gaseous disk 
becomes neutral (independently of the passive-disk outcome advocated in Menou et al. 2001). 



3.3. Observational Constraints on Disk Sizes 

The presence of disks around pulsars cannot be hidden: even a relatively compact disk will 
produce some optical emission, as radiation at these wavelengths is produced at relatively small 
radii. 

Optical emission is a common feature of pulsars: since the first detection of the Crab in 1969, 
several more pulsars have been observed in the optical bandpass, and some of them have been 
seen to pulsate. In the cases where no pulsations are seen, the optical radiation is believed to 
be mainly of thermal origin, while pulsed optical emission is generally considered to be either a 
mixture of thermal and non-thermal (magnetospheric) radiation, or purely magnetospheric. Pacini 
(1971) and Pacini &: Salvati (1987) proposed that the high energy emission comes from relativistic 
electrons radiating via synchrotron processes in the outer regions of the magnetosphere. 

The Crab, Vela, and PSR B0540-69 pulsars belong to the above class, all possessing a 
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pulsating optical counterpart. On the other hand, no pulsations have been detected for the 
candidate optical counterpart of the source PSR B1509-58 proposed^ by Caraveo et al. (1994). 

In the following, we estimate, for each of the sources that we consider (except for PSR 
J1119-6127, for which no optical limits are available so far), the maximum size allowed for a 
disk to be compatible with the observed optical limits. For the three sources for which the 
optical emission is pulsed, we consider as an upper limit the lowest value of the emission (i.e. the 
minimum of the light curve), while for PSR B1509-58 we consider the detected value, assuming 
that it is indeed associated with the pulsar. 

The inner radius of the disk is assumed to coincide with the light cylinder radius, Ric- The 
disk model that we adopt is similar to that described by Perna, Hernquist & Narayan (2000). We 
include both the flux due to viscous dissipation and that due to reradiation. The irradiating X-ray 
luminosity is assumed to be isotropic, because, due to light bending in the vicinity of the star, 
beaming effects are smoothed out (Perna & Hernquist 2000). 

Using the X-ray luminosities and periods from observations (see Becker & Trumper 1997 for 
a review), we find (assuming a disk inclination angle of 60°) that, for a mass inflow rate ~ 10^^ g 
s~^ a disk of radial extent ^ 6 x 10^ cm is compatible with the optical limits for the Crab pulsar. 
A large disk (^ lO^^ ^m) is allowed for both PSR B0540-69 and PSR B1509-58 at the mass inflow 
rate of ~ 3 x 10^^ and ~ 10^^ g s"^ respectively, while for Vela only a very small disk (^ 10^ cm) at 
an high inclination (^ 89°) would be compatible with the observational constraint^. The allowed 
size becomes larger as M decreases. Again, we note that, with the exception of the Vela pulsar, 
these disk sizes appear consistent with the values expected at late stages (before the fallback disk 
becomes neutral) in the evolutionary scenario of Menou et al. (2001). 

Support for this work was provided by NASA through Chandra Fellowship grant PF9-10006 
awarded by the Smithsonian Astrophysical Observatory for NASA under contract NAS8-39073. 
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